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Abstract 



Cn , The method of multiscale analysis is constructed for dicrete systems 

of evolution equations for which the problem is that of the far behavior of 
^O ' an input boundary datum. Discrete slow space variables are introduced 

in a general setting and the related finite differences are constructed. The 
ZTi ' method is applied to a series of representative examples: the Toda lattice, 

(«-»* \ the nonlinear Klein-Gordon chain, the Takeno system and a discrete ver- 

qv . sion of the Benjamin-Bona-Mahoney equation. Among the resulting limit 

Q\ ' models we find a discrete nonlinear Schrodinger equation (with reversed 

space-time), a 3- wave resonant interaction system and a discrete modified 

Volterra model. 



O '. 1 Introduction 

C*3 



The study of nonlinear dispersive waves, historically originating in water wave 

problems, has received considerable attention as a consequence of two funda- 

jJJ ■ mental discoveries. The first one is the concept of complete integrability of 

partial differential equations discovered for the Korteveg-de Vries equation jy 
and soon extended to the nonlinear Schrodinger |2| and sine-Gordon J3(] equa- 
tions (for a review see j|, |5|, || ) . The second one is the modulational instability 
of wave trains [j/J which is a mechanism for the creation of localized nonlinear 
solitary waves, the solitons ||. 

These two discoveries have a common origin, the reductive perturbation 
method (or multiscale analysis) H allowing to deduce simplified equations from 
a basic model without loosing its characteristic features. The method consists 
essentially in an asymptotic analysis of a perturbation series, based on the exis- 
tence of different scales. More specifically, the method generates a hierarchy of 
(small) scales for the space and time variations of the envelopes of a fundamental 
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(linear) plane wave and all the overtones. The scale is moreover directly related 
to the (small) amplitude of the wave itself. The scaling of variables is performed 
via a Taylor expansion the frequency u>(k) in powers of a small deviation of the 
wave number fc = feo + eK - This deviation from the linear dispersion relation is 
of course generated by the nonlinearity. 

The success of the method relies mainly on the nice property that the result- 
ing reduced models are simple, representative and often integrable. Simple here 
means actually simpler than the master equations and allowing for usefull in- 
formations. Representative means that they illustrate effectively real processes. 
This property relies on the self consistency of the perturbation series which 
treats all overtones and avoids secularities |1C| , |llf| . Integrable means that they 
carry an infinite set of conserved quantities, have a bi-hamiltonian formulation, 
are solvable (in some sense), etc... Finally, as emphasized in |1J], there ex- 
ists a general property of the reductive perturbation approach which allows to 
understand, in a qualitative way, why the reduced systems are often integrable. 

The situation is quite different in the case of nonlinear lattices (continuous 
time and discrete space) for which a reductive perturbative method does not 
exist which would produce reduced discrete systems. There are actually three 
different approaches to multiscalc analysis for a discrete evolution. The first one 
is obviously to go to the continuous limit right in the starting system, for which 
discreteness effect are wiped out. The second one is the semi-dicrete appoach 
which consists in having a discrete carrier wave modulated by a continuous 
envelope |13|. In that case some discreteness aspects are preserved, in particular, 
the resulting modulational instability may depend on the carrier frequency. The 
third one stems from adiabatic approximation, but the approach requires to use 
the rotating wave approximation to artificially eliminate the overtones. The 
price to pay is that the predictions, e.g. the modulational instability, are not 
trustworthy for large times [ H . 

We propose here a set of tools which allow to perform multiscale analysis on 
a discrete evolution equation when the problem is that of the propagation of a 
signal sent at one end of a nonlinear lattice. These tools rely on the definition of 
a large grid scale via the comparison of the magnitude of the related difference 
operator, and on the expansion of the wave number in powers of frequency 
variations due to nonlinearity. 

The method will be illustrated in a series of examples (Toda lattice, anhar- 
monic Klein-Gordon chain, BBMP equation, ...) for which the reduced model 
for the slowly varying envelope ijj(n, t) results in the following evolution (over 
dot stands for partial time-derivative) 

-i(3 [ip n+ i - Ipn-l] + otij) n - 7|^>„|Vn = . (1.1) 

The coefficients a and [3 and 7 depend on the starting model equation and on 
the frequency of the carrier wave. The continuous version is well known to apply 
to pulse propagation in a nonlinear Kerr medium (optical fiber) |lq] and has been 
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obtained also in the context of Rayleigh- Taylor instability and electron-beam 
plasma Eg ] and refered there as the unstable nonlinear Schrodinger equation. 

Another example deals with the Takcno model of exciton-phonon coupling in 
diatomic chains |L7| for which we investigate the three wave resonant interaction. 
The resulting discrete system reads 



(1.2) 



where X is the envelope of the phonon wave (carrier frequency f2, group velocity 
v), aj are the envelopes of the two components of the exciton wave (carrier 
frequencies u- r group velocities Uj), for the Brillouin selection rule V\ — im = H. 

We illustrate also the method by investigating the low frequency discrete 
limit of the Toda lattice [pi , and obtain 

U n +1 - Un-l + U n U n = , (1.3) 

which can be viewed as the Volterra model JlSJ ] with space and time exchanged. 

In this paper, the emphasis will be put mainly on the method and on the 
above mentionned simple examples, but not on the physical implications that 
one can infer from the reduced systems. Indeed, such studies essentially depend 
on the physical context and thus require specific attention. After a short state- 
ment of the problem in the next section, the sec. |3| is devoted to the definition 
of the main tools. The sec. [| deals with the slowly varying envelope limit of the 
Toda lattice, and sec. g with nonlinear Klein-Gordon (or sine-Gordon) chain. 
The sec. o deals with a new nonlinear evolution constructed as a discrete version 
of the Benjamin-Bona-Mahoney-Peregrini equation, and the sec. [7] is devoted 
to the study of three wave resonant interaction in the Takeno model. 



2 Discrete waves 

2.1 Boundary value problem 

The pupose is the study of a nonlinear dispersive chain with dispersion relation 
fl(K). The physical problem we are concerned with is the following: the first 
particle of the chain (say n — 0) is given an oscillation (or is submitted to an 
external force) at frequency O. Would the chain be linear that this oscillation 
would propagate without distorsion as the plane wave exp[iftt + Knd)], with d 
being the lattice spacing. 



But the nonlinearity induces some deviations from the value fi, namely, the 
wave propagates with actual frequency u> and wave number k that we define as 

LU = Q + ev, (2.1) 

and by Taylor expansion 

k = K + e - v + e 2 7 j/ 2 , (2.2) 

c 

for the following definitions (note that c is the group velocity at K) 



1 _ dk 
c duj 



- — 
duj 2 



(2.3) 



For notation simplicity, we shall assume everywhere that 7=1, which does not 
reduce the generality of our task. 

A wave packet in the linear case is then given by the Fourier transform 

Un(t) = I ' du u{uj)S" Jt+knd \ 



where u(u) has support inside the allowed band frequency [wo,W£,], which gives 
here 

u n (t) = e S nt + Knd ) J dv u(v)e lve{t+nd/c) e ™ 2f2nd 

or else 

r A(n, t) = S nt+Knd ^ 
Un(t) = A(n,t)ip(£ n ,T n ), I (2.4) 

{ iP{Zn,T n ) = efdvu(v)e^+" 2 ^ 

by means of the following change of independent variables 

r n = e(t + nd/c), Cn = e 2 n. (2.5) 



2.2 Discrete scaling 

In order to keep discreteness in the space variable for the envelope ip(£, n , T n ), we 
fix the small parameter as 

e 2 = 1/JV , (2.6) 

and, for any given n, we shall consider only the set of points {• • • , n — N, n, n + 
TV, • • • } of a large grid indexed by the slow variable m, that is 

• • • , (n - TV) -> (m - 1) , n-> m , (n + N) -> (m + 1) , • • • (2.7) 



To simplify the notations, we shall be using everywhere 

^(£/.7"i) = i>i > ^(C/» T n) = ^j , (2-8) 

for a given n and all j (note that ip n = ip n ). Hence we are interested in expressing 
everything in terms of ip m = tp(m, r) defined as 

tpn-N = lpm-1 , 1pn=1pm , 1pn+N = VWl- (2-9) 

The problem is now to express the various difference operators occuring in non- 
linear evolutions for the product A(n,t)ip m in terms of difference operators for 
ipm- 



2.3 Initial value problem 

The traditionnal approach to multiscaling for continuous media originates from 
water wave theory for which the physical problem is ususally that of the evolu- 
tion of an intial disturbance (e.g. of the surface). In that case the observer has 
to follow the deformation at the (linear) group velocity. 

This operation corresponds to making, in the general Fourier transform solu- 
tion, the expansion of u>(k) around small deviations of k from the linear disper- 
sion law. The resulting change of variables (in the discrete case for comparison) 
would then read 

T„ = e 2 t, £„ = e{nd + ct) , (2.10) 

which indeed corresponds to a translation at the group velocity c, that is in the 
co-moving frame. 

In that case, the change of variable definitely breaks the discreteness of the 
space variable. Hence an initial value problem for a discrete lattice in continuous 
time cannot be treated within a fully dicrete multiscale analysis. 

However, when the phenomenon to observe results from a boundary input 
datum, the observer stands at some given point of the lattice and compares 
the received signal with the input signal. This has to be done in the retarded 
time, which corresponds to the change of variables ( |2.5| ), and obviously does not 
destroy the discrete character of the variables. 



3 Difference operators 

3.1 General consideration 

Before going to the technical points, it is quite usefull to recall elementary facts 
concerning the relation between continuous derivation and discrete differences. 
Let (p n denote a function of the discrete variable n which varies slowly from one 



site to the other. In that case we can represent this function by the function 
f(x) = 4>n of the continuous variable x — nd. The constant d is a small arbitrary 
real number which needs not to be related to the dimension of the grid spacing 
(this would be necessary if f(x) would be the datum and cj> n its representation). 
By Taylor expansion we readily obtain 

d ^(x) = _^ + 2 

ox 

d djp_ = 1 _ + Q 3 

02: 2 

d2 ^-Qr = ^+i~ 2 ^ + K-i+o{d A ) . (3.3) 

Consequently, if the precision of the first derivative has to exceed at least the 



value of the second derivative, it has to be defined with (3.2). This is not the 



only reason for considering centered differences. The first order wave equation 

dt4>n - V -[4>n+l - 4>n-l] = , 

possess a good dispersion relation (w = v sin fc) with real valued solution. Such 



would not be the case with the definition (3.1) 



Going now to the stretched space variable £ = ex, we readily get 
+"*-+» ±e d -8T +e 2* -d?- ±e 6* -8p- + '" 



which, under discretization procedure according with rules (3.2) (3.3), lead to 
the following qualitative indications for the change of variable on the first and 
second derivatives 

<£n+i - K-i = e 2 [^m+i ~ 4>m-i] + 0{e 6 ) (3.4) 

4>n+l - 20„ + <f) n -l = e 4 [0m+l - 2<f>m + 0m- 1] + C(e 8 ). (3.5) 



3.2 Construction of the stretched grid 

For e 2 = 1/N, the above expressions serve as a guide to the correct rules for the 
chan ge of coo rdinates. They can be understood simply within the definitions 
(U) and (|3) b y writing 



4>n+N — <f>n-N = {.4>n+N — ^n+N-l) + {4'n+N-2 — <frn+N-4) 

+ • • ■ + (0n-JV+4 — <f>n-N+2) + {4>n-N+2 — ^ji-n) 
~ N[<j> n+ i - <j> n -l] 

where N has to be chosen odd such as to find in the above expansion the very 
term <fi n +i — <j>n— l- Then the meaning of this grid change is that, on the interval 



[n — N, n + N] , the variations of <f>„ are almost equal, or else that <j> n is of almost 
constant slope (with a precision of 1/N 2 ). The purpose now is to construct 
rigorously such simple rules, and in particular to determine the conditions of 



slow variation under which expressions like (3.4) and (3.5) do hold 



To that end we define the derivatives in the original variable n as (the C| 
arc the binomial coefficients 



e\(k-ey.> 



V0„ = <f>, 
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vV« = J2(-Y ci 



■>n+k-2£ , 



(3.6) 



and the derivatives in the new varaible m defined in (2.7) as 



A 



AK 



\-N 
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n+N(k-2£) 



(3.7) 



£=0 



By expressing the value of a function <f> n at any other point n + N in terms 
of the differences V fe (/>„ up to k = N, we prove Q in the appendix that, for 
any odd N = 2q + 1 



±n(, 



^■NVn 



£a^V 2 ^Vn, 



2q+l 



E < v2 ^» . 



,_ (2g+l)(g + l)l 
9 («-/)!(2€+l)! ' 

< 2(2g+l)(2 g + £)! 
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(2g + l-£)!(2*)! 



(3.8) 



(3.9) 



The above expressions are identities which replace the analogue continuous 
change of variable d x <f> = e 2 d^(f>. Note that the first coeffcients are a® = N and 
7,i = TV 2 , so as for the p-th order A at difference whose first term is N p V p 4> n . 



3.3 Slow variables 



We make now the hypothesis of slow variation of the function <j> n as 

|V fc+1 0n| = e 2 |V fe 0„|+O(e 4 ), (3.10) 

for some norm. Such an hypothesis is never explicit in the continuous case where 
it is usually admitted implicitely that if d x (f> — e 2 d^cj), then \d 2 (f>\ ~ e 2 |c*£ C <^ > | - But 
actually it goes the other way round: first one must assume that the function 
4>(x) is of slow variation in x, i.e. 



d k+1 4> 
dx 


~e 2 


d k (j) 
dx 



and second it is this very measure e 2 of the velocity of variation of <j>{x) that 
allows to make the change of variable £ = e 2 x. 



It is necessary now to check that the hypothesis of slow variation is compati- 



ble with the expression (3.8) and (3.9) as soon as N is chosen equal (nearest odd 
integer) to e~ 2 . By a recursive use of (3.10) we obtain the following majorations 
at large N = 2q + l: 



^|Am 
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|v^„|£ 



o , 
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■]y2\ A N<t>n\ 
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iv 2 ^„i J2 



(2q + 1) 2 ^+1 
>1 



t (2<7+l) 2 * ' 
We compute the sums of the above series as N = 2q + 1 — v oo : 



limV 



, -*> — (2q+l)"- 



t = E 



e=o 



2q+l 

lim y 



, -» — (2g + l)2^ 



E 



-— - = 2sinh(-) = 1.042 

2 2 ^+l)! V 



= 2cosh(l) -2 = 1.086 



(21)! 



(3.11) 



(3.12) 



Hen ce the slow variation assumption are compatible with the series ( |3.8|) and 
(3.E). The convergence is actually quite fast as, e.g. in the order 1 difference, 
the first correction coefficient is 1/24 to be compared to 2sinh(l/2) — 1 (they 
differ by 0.5 x 10" 3 ). 



Now, with the assumption ( 3.10 ), the identities (3.8) and (3.9) together with 
1/JV = e 2 , N = 2q+ 1, lead to 



V<£ n = - A N n + O{l/N 3 ) , V 



- ■ = — AUn+0(l/N i ), (3.13) 



which are the general rules for the change of variable n — > m in the first and 
second differences. These rules do concur with (EO) and tpM. 



3.4 General expressions 



Remembering the notation (2.8), by Taylor expansion in the continuous variable 
r„, we have 

ijj n+1 = ip n+1 + e-d r ijj n+1 + -(e-) 2 (9 2 V>n+i + -{t-fdli> n+ i , 
c 2 c be 

i'n-i = ipn-i - e-d T ip n -i + -(e-) 2 d?Vn-i - -(e-) 3 9 3 V«-i • 
c 2 c be 

Then we use the identities (such general expressions are given in the appendix) 

lpn+1 - Ipn = ^bPn+l - 1pn-l] + ^ bPn+1 - 2 ^n + i>n-l] , 
Ipn - 1pn-l = x[V^+l - 1pn-l] ~ ^ [V>n+1 ~ 2lp n + lp n -l] , 



to replace hereabove ip n +i and ipn—i m terms of ip n . Finally by means of (3.13) 
we arrive at the following general expressions which constitute the basic tool 
allowing to compute the first and second derivatives: 

ip n+1 =4> m + e-d T ip m + -e 2 (-) 2 dlip m + — -[ip m +i ~ ipm-i] 
c Ac iv A 

1 e 2 d 

+ 2jy2 WWl - 2 ^m + tpm-l] + —(-) d r [lp m+1 - 1p m -l] 

+ ^(-/d% n + 0(e 5 ), (3.14) 

$n-i =ipm- e-d T ip m + -e 2 (-) 2 d 2 tp m - —-[ip m+1 - ip m -i] 

C Ac iV A 

N CA DC 

1 e 2 d 

+ -^p. bPm+l ~ tym + Ipm-l] ~ Tjy(~) d r[i>m+l ~ tpm-l] 

+ ^(- c ) 4 d%n + +0(e 5 ). (3.15) 



Note that pushing the general formula ( 3.13 ) to higher orders, we can reach 
any order of precision in such expansions. 



3.5 Discrete derivative of a function 



By straightforward calculations, using the above basic formulae ( [3. 14 ) and 
( |3.15 ), the first derivative follows 



4>n+l ~ $n-l = 2e-5 T ^ m + —[lp,n+l ~ 1pm-l] + 0(e 3 ) (3.16) 

c iv 

together with the second derivative 

$ n+ i-2$„ + $ n _! = e 2 (-) 2 d 2 ^ m + — - d T [i> m +i - ip m -i} 

c IN C 

+ T^m+i - Wm + */> m _x] + L.(-)*%il> m + 0{e 5 ). (3.17) 

iv z 12 c 

It is worth remarking that, by continuation of the above formulae stopped 
at e 3 for the first derivative and at order e 4 for the second we obtain 

d _ 2 & e & n( 3\ 

dx d£ c dr 

q2 e 3 d d e 2 d 2 

dx 2 c dr <9£ c 2 dr 2 



which are the very relations corresponding to the change of variable 

t = e(t + x/c), £ = e x , 



continuous analogue of (2.5) 



3.6 Discrete derivative of a product 

We have now to obtain the analogous relations for the product 
u n (t) = A(n,t)ip(£n,T n ) = A n ij) n 



appearing in the definition (2.4). The quantity w„+i — u n _i is factorized as 



A n+1 1p n +l - An-ltpn-l = 

-[A n +i - A n _i][^„ + i + i) n -i] + -[A n+ i + A„_i][Vj„+i - ij>n-l\ 



and using the basic formula ( |3.14 ) and ( 3.15 ), it follows, at order e/N or e 



u n+1 - u n -x =[A n+1 - A n -i]ip m + e-[A n+1 + An-i]d T ip„ 

c 

— -[An+x + A n -i][lp m+ i - 1p m -l] + 

e 2 {-) 2 ]-[A n+1 - A n ^]d 2 T ^ m . 
c 2 

The second derivative is factorized as follows 



(3.18) 



U n+ i - 2u n + U n -i = ~[A n+1 - An-i][lp n +i - V>„-l] 

+ ~[A n+1 + An-l][i)n+l - 2"0n + 4>n-l] + 4>n[A„,+l - 2A n + An-l] 



which readily gives from (3.14) and (3.15) 



ln+1 - 2u„ + U n -1 =[A n +\ - 2A n + A n -l] "0m 

+ e[A n+ i - A n -i]-d T il> m 
c 

+ ^[4 +1 + A, l _ 1 ](-) 2 ^ ra 
2 c 

+ -rr-[A„ + i - An_i][lp m+ i - 1p m -l] 



(3.19) 



at order e/N or e 3 . The fromulae ( |3.18 ) and ( 3.19 ) constitute our basic tool for 
deriving reduced models in the next sections. 
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4 The Toda lattice 

The Toda chain is defined by the equation p8| 

x n = e x "+ 1 - x --e x "- x - 1 , (4.1) 

which can also be written 
dB{n,t) 



dt 
dV(n,t) 



[1 + B{n, t)]{V{n, t) - V(n - 1, *)), (4.2) 

B(n + l,t)-B(n,t) (4.3) 



dt 

for the following definitions 

1 + B(n,t) = e x "- X "- 1 , V(n,t)=i„. (4.4) 

The dispersion relation for the linearized version of the above equation 

ft 2 = 4 sin 2 — (4.5) 

is precisely that of the chain of coupled masses. 

Following section |2j we expand oj = il + tv and consider here the two cases 
Q = and f2 ^ 0. The case J7 = can be viewed as the low frequency limit of 
the Toda chain, which is actually the standard long wave limit. The case $7 ^ 
corresponds to the slowly varying envelope approximation of solutions of the 
Toda lattice. 

4.1 Slowly varying envelope approximation 



We first seek a solution of equations ( |4.2| ) and ( |4.3| ) in the form of a Fourier 
expansion in harmonics of the fundamental A(n,t) — expi(£lt + Knd) where 
the Fourier components are developped in a Taylor series in power of the small 
parameter e mesuring the amplitude of the initial wave 

l=p p=oo 

B(n,t) = J2 £ ^i l) (Cn,r n )A l (n,t), (4.6) 

l—~p p—1 

l—p p=oo 

V(n,t) = J2 E e P 4' ) (&.'-»)A I (n,t). (4.7) 

i——p p—i 

Note that the above series includes all overtones A ; (n,i) = expil(£lt + Knd) 
up to order p. These are generated by the nonlinear terms which explains that 
the corresponding coefficients are of maximum order e p . Here we have the real- 
valuedness conditions 

^-° = (^°r, 4~ l) = (4°r (4.8) 
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the asterisk denoting complex conjugations. The slow variables r n and £ n are 
introduced via 



nd 2 

c 



(4.9) 



where the velocity c will be determined later as a solvability condition of equa- 
tions (O) and (fO). 



By substitution of (|4J) tyLty into (|4j) and ^4J) and using ( |3 . 1 8| ) (together 
with dt = td T ) at order two in e we obtain 



1= p p— oc 



a 



E E ^TT^ ("», r) + ifity« (m, r)} A' (n, i) 



l = —p p=l 



9r 



(4.10) 



i 1+ E E^W^W)} 

i——p p—i 

l—p p—OO 

X <E E eP [4°KT)(^(n,i)-^(n-l,t) 

l——p p—1 

+ y4"K^ , l e 2 d ^p'W) 

c or 2 c OT 2 

+ 6 2 (^)(m + l,r)-^)(m-l,r))^(n-l,t) + ...]} (4.11) 



l—p p—OO 



E E <?{*»-*? {m > T ) + m $ ) ( m ' r » A ' ( n ' *) 



i——p p—i 



dr 



l=v 



£ e p {4 l) (m,T)(A\n + l,t)-A l (n,t)) 



/= 



-,(0/ 



C OT 

^(y 2 ^ ( ;- T) A'(, + M) 



2 c 



cV 2 



+ e 2 (^(m + 1, r) - ^W(m - 1, r))lA'(n + 1, i) + ...} (4.12) 



We can procee d now to co llect and solve different orders of e p and harmonics 
I, order (p,l), in ( 4.10 ) and ( 4.12 ). Note that it is e nough to consider I > as 
negative values follow from the reality condition (4.8). In the leading order (1,1) 
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we have 

1=1 1=1 

Y^ Kilt/tip (m,T)A\n,t)- J2 <t>i\m,T){A l {n,t) ~ A l (n - l,i)) = (4.13) 
i=-i i=-i 

i=i i=i 

J2iM<j) { p(m 7 T)A l (n,t)- ^ ip[ l) (m,T)(A l (n+ l,t) - A l (n,t)) = (4.14) 

This is a linear homogeneous system for ip^ (to, r) and 0^ (to, r) polynomial in 
A. Hence each coefficient has to vanish separately. For I = the system gives 
trivial equations. For I = 1 the determinant of the system for i/j\ '(jn,T) and 
4>\ (to, t) is zero if Q(K) verifies the dispersion relation 

Kd 

n(K) = 2sm— (4.15) 



Under this condition we seek the general non-trivial solution of equations (4.13) 
and ( |4.14| ) 



ij)y (to, t) = ar)(m, r) , 0j ' (to, t) = 77(771, r) (4-16) 

with a given by 

,—iKd. 
a = exp(^— ) (4.17) 



At order (2,/) we have the system 
£ {ifirfKr) - ^°(m,r)(l - o 2, )}A'(n,*)] 

l=-2 

+ E { ^ m ' T) - (*)M^ * }i4 « (n , t)] 
J=-l 
1=1 i=i 



{J24Hrn,r)A l (n,t)x £ <^W)(1 - a 2l )A l (n,t)} = (4.18) 



;=-i ;=-i 



J{i(lJ4 I) (m l T)-4 ,) (m l T)(a- J '-l)}A , (n,t)l 



Z=-2 



+ g { d£ V, r) _ ( d } g^ V, r) (o -2« M , (n> t)] = (419) 

Z=-l T C 

For / = we obtain an homogeneous system with non-zero determinant, conse- 
quently only the trivial solution exists, so that 

V>i 0) (m, r) = , <^ 0) (to, r) = (4.20) 
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For I = 1 we have an inhomogeneous linear system for ijy^ ( to j t ) an d 4>2 ( to j T )- 
The determinant of the associated homogeneous system is zero owing to the 



dispersion relation ( 4.15 ). Therefore, the system will have a solution if the 



Fredholm solvability condition is satisfied, that is if 



dVL(K) Kd ,,„,, 

c=^i = rfcos- (4.21) 

which determines c as the group velocity. 
Under this solvability condition we get 

,(■11/ -, 1-/ n a , d .driim.T) ,, „„. 

</4 V> r ) = aS{m,T) + —{a- - 1) '^ J (4.22) 

$\m,T)=8(m,T) (4.23) 

where S(m, r) is an arbitrary function. Furthermore for I = 2 we obtain 

1p2 ( TO : r ) an d </>2 ( m i T ) 

2(1 2^ — sin Kd) 



(2), ^ siri(Kd) 2/ 

2(f2 2 -sin 2 ifd) r 



^ J (m,r) = — — 4 2 l^ 2 (m,r) (4.25) 



a 2i 



The next order (3,1) gives the system 

1=3 

J2 {im4\m,T) - 4V,t)(1 - a 2l )}A l (n,t) 

l=~3 

^ f dTp [ 2 l \m,T) _ d d^Hjn^r) 2l ,l,A 2 9^iW) 
+ ^ { Or V dr a + W 9r 2 

l = -2 

- [^\m + l,r) - 0l°(m- l,r)]^}A'(n,t) 

Z=l Z=2 

- ^ vi I) (w*,T)i4 , (n,t)x XI ^K^l-a^M) 

Z=-l Z=-2 

1=2 1 = 1 

- Y, 4 l) (rn,r)A l (n,t)x £ ^\m,T)(l - a 2l )A l (n,t) 

l = -2 l=-l 

- £ 4 l \m,r)A l (n,t) x ]T (-) d<l>1 ^ T a 2l A l (n,t) = (4.26) 
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;=3 



£ {irW^V, t) ~ 1$V> r)(a- 21 - l)}A l (n, t) 



2=-3 



g { MV,r) _ ( y4'W) fl - a}i , (M) 

2=-2 



g ( i )( g )a gfelO oMi4 , (w>t) 



k 2 ; V 9r 2 

2=1 



1=1 22 

2[^' ) (m+l,r)-V? ) (m-l,r)]^ , (n,t) = (4.27) 



The order (3,0) allows us to determine ip2 \ m i T ) an d </>2 (to,t) via an inho- 
mogeneous system of equations. They read 

4° ) (m,r) = -^|r 7 (m,r)| 2 (4.28) 

^'K^^llK^f (4-29) 



The next order (3, 1) is a tedious one which results in the nonlinear evolution 
for T](m,T). It is an inhomogeneous linear system of equations for ip^ (to, t) 

and cf> 3 (to, t) of determinant zero. Therefore, we will have a solution if the 
Fredholm solvability condition is satisfied. This condition gives a nonlinear 
evolution of 77(771, r) in which the term in 6(m, t) coming from 1P2 ( to j T ) an d 
4>2 (ttt-, t) is self-eliminated. The equation for 77(777, t) = ?7 m finally reads 

f) 71 

-7/3[?7m+l - Vm-l] + a ^TT ~ 2 \Vm\ 2 Vm = , (4.30) 



with 



Kd 
a = tan 2 , [3 = sin Kd . (4.31) 



4.2 Low frequency limit 

We are interested now in the propagation of a low frequency wave and we define 
the following quantities 



duj 
UJ = 6 "> -dk 



= v g = d, (4.32) 

fc=0 



and then we develop the function k(u>) in Taylor series in e. We obtain 

1 , , 1 d 3 k 



fc = _ w + 7e 3 l/ 3 ; 7= _ ^ 
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(4.33) 

fc=0 



From the above expansion, the new variables follow 



T„ 



■„ = e(t + n) , Cn = e n, (4-34) 



and consequently we shall consider the new grid with the definition 

„ 1 

e =— , £„ -> m , C«+jv -> m + 1. (4.35) 

The method and formulae developped before applies identically except that 
now we have \/N — e 3 instead of e 2 . In particular, the only formula to use is 



that of the second derivative (3.17) 
Now we set 



B n (t) = e 2 b(U,T n ), (4.36) 



and compute the limit expression of the Toda equation (4.4). We obtain 



^ ( tT^T ^ = £ 4 d TT b(m) + e 6 d T [b(m + 1) - b(m - 1)] + G(e 8 ), 
\1 + e z b{m) J 

which finally reduces to 

b(m + 1) - b(m - 1) + b(m)d T b(m) = 0. (4.37) 

This equation is a Volterra-like equation where one would have exchanged space 
and time, as expected. 

5 Nonlinear Klein- Gordon chains 

The modulation of the solutions of the Toda chain has been shown to obey the 



NLS-like equation (1.1 ) with particular values of the coefficients. It is of interest 
to compare this situation resulting from an integrable model to the one resulting 
from a non- integrable starting equation. Such is the case for the nonlinear 
Klein- Gordon chain 

U n - Uj'l (U n +1 - 2u„ + U„-l) + LU^Un + Tu^ = , (5.1) 

or the sine-Gordon chain 

ii„ - lo\ (u n+1 ~ 2u n + u n -i) + ^o sin u n — • (5.2) 

Both cases, in the perturbation scheme, are equivalent for T = — 0Jq/6, and their 
dispersion relation is 

tt 2 =c^+4u; 2 sin 2 (^), (5.3) 
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5.1 Evolution of the envelope 



We start with the evolution ( |5.l| ), and seek the evolution of ipm = V'i with the 
tools developped in sec. H when 

oo i=p 

u(n,t) = Y. eP E e m ^*V P {m,T) (5.4) 

p— 1 £——p 

with conditions on ^ which ensure reality of w n : 

V^ = ^ , ^ 0) € M. (5.5) 

Here 9(n,t) — fit + Knd, the slow variables are those defined in sec. (or 



those used for the Toda chain, namely (4.9)), and we stop everything at order 
e 3 = e/N. 

The coefficients of the constant term give at order e : 

wjty? =* < = , 
atrder e 2 : 

and at order e : 

«) TT -c 2 (^) 2 «) TT + r[(^) 3 + 6|^T^] + ^3° => ^3 = o. 

The coefficients of e ie at order e give 

^J[-fi 2 - w 2 (e iKd - 2 + e- jA ' d ) + lo 2 } = , 

which agrees with the dispersion relation ( |5.3| ). Then at order e 2 we obtain 

i\}\ [-CI 2 - u 2 {e lKd ~ 2 + e~ iKd ) + u 2 ] + 

(iPl) T [2ifl - u 2 -{e lKd - e- lKd )} = 0. 
" ~ c 

The first line here above vanishes due to the dispersion relation while the second 
one vanishes too thanks to 

n = uj 2 - sin Kd. (5.6) 

c 



which is indeed verified as, by the definition ( |2.3| 

dfl 1 



c = fk = h uj2id siaKd - ( 5J ) 
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Finally the order e leads to 



%l>\ [-VL 2 - uj\ [e lKd - 2 + e- lKd ) + J 2 ] + 

(^i) T [2»n - ^-(e lKd - e - jKd )] + 
c 

(^) Tr [l- w ^(^) 2 (e^ + e-^)j 
2 c 

- ^~(e^ d - e- iKd M{m + 1) - <#(m - 1)] + 3T|^| Vi- 

The first two lines hereabove vanish identically and we are left with the equation 
for ip = tpl, 

-i(3 [Vwi - VVn-i] + aip TT {m) + 3T\Tp m \ 2 i/: m = 0, (5.8) 

with the definitions 

(3 = uj\ smKd , a = 1 - ojf(-) 2 cosKd. (5.9) 

c 

The main difference between the parameter values hereabove and those for 



the (integrable) Toda chain ( 4.31 ) is that here a changes sign for some value 
of K, which is of fundamental importance for the stability properties of the 
modulation. But this problem goes beyond the scope of this paper and will be 
considered elsewhere. 



5.2 Continuous limit 

In order to check the consistency of our method, it is instructive to exam- 
ine the continuous limit of the nonlinear Klein-Gordon chain to verify that its 
(continuous) multiscal e an alysis gives rise to an equation which is precisely the 



continuous version of (5.S) 



Defining the continuous variable x — nd and the velocity v = tvid, the 



continuous limit of (5.1) reads 



u u — v u xx + UJ u + Tu' = 0. (5.10) 

By seeking a solution as a Fourier integral and expanding k in powers of e for 
lo = fi + ev around K, we are led as previously to the expansion 

u( Xl t) = Y, eP E ^'''Wp&r) , e = nt + Kx, (5.ii) 
p=i i=—p 

and the change of variable (c is the group velocity at frequency f2) 

Z = e 2 x, r = e (*+-). (5.12) 

c 
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with the reality conditions ipz, = ipp 



Inserting everything in the evolution equation (5.10) we get as the coefficients 



of, 



0=e[-n 2 + v 2 K 2 +uj 2 ]+e 2 



2itiip T - 2iK—ip T 



f T - 2iKv 2 ^ - — Vvr + 3TA\^\ 2 1JJ 



0(e 4 ) 



The order e cancels as soon as we select the linear dispersion relation 



n 2 = 



v 2 K\ 



(5.13) 



The order e 2 then identically vanishes as indeed, from the definition of the group 
vclocit, we readily get 



K v 2 

c = v 2 — => ifl iK = 0. 

it c 



(5.14) 



Finally, the order e 3 furnishes 



-2iKv 2 ^ + [1 - -]i/> TT + 3r|vf V> = o. 



(5.15) 



One interesting consequence here is that the coefficient of the second derivative 
hereabove is from (|5.14) 



i_f S! i_i£ = _J!L < o 

9 9 tv-9 9 r^9 ' 



v 2 K 2 



l K 2 



(5.16) 



which never changes sign, contrarily to the discrete case. 

The continuous equ ation ( 5.10J ) being the continuous limit of the discrete 
Klein-Gordon equation (5.1), the consistency c heck consists now in the verifica- 
tion that the continuous limit of our equation (5.S) is precisely the above NLS 
equation ( 5.15 ). This readily follows from the limits as Kd — > : 



r[^m+i - Tpm-l 



-> d ^ 
u>id — f v 



a= 1 



(3 = lo\ sin Kd 
t. — cos Kd 



,K 



6 Discrete Benjamin-Bona-Mahoney-Peregrini 
equation 
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6.1 Construction of the model 

From the quite well known Boussinesq equations p0| , further asymptotic lim- 
its and restriction to unidirectional propagation are possible, leading to re- 
duced models among which we find the Korteweg-de Vries equation (KdV) and 
the following Benjamin-Bona-Mahoney-Peregrini (BBMP) equation for the field 

u(x, i): 

u t + u x - u xxt + uu x = 0. (6.1) 

The discrete multiscale tool, unlike in the continuous case, furnishes expres- 
sions for the differences as infinite power series in the small parameter. It is 
always more convenient to stick with the first few orders in the expansion, like 
in the basic expression ( 3.16| ). Consequently we need to get rid of the third 



order derivative in the above evolution by writing instead ( |6.1[ ) as a system of 
equations for u(x,t) and two auxilair fields p(x,t) and v(x,t) 

u x = P, 

Pt = v, (6.2) 

v x - u t =p(l + u). 

The discretization of a continuous model is always non unique. For instance, 
the discrete analogues of KdV are the Toda lattice or the Langmuir equation, 
both leading to KdV in the continuum limit. Among different possible choices, 



weintroduce here the following discrete analogue of (6.2): 

-[u(n+ l,t)-u(n- l,i)] =p(n,t), 

dp(n, t) 

— Ol — =v(n,t), (6.3) 

l[v(n+l,t) - v(n-l,t)} - d^*) = p ( n ,t)[l + u(n,t)]. 



6.2 Multiscale analysis 



As in the Toda case we seek a solution of equations (6.3) in the form of a Fourier 
expansion in harmonics of the fundamental A(n,t) — expi(£l(K)t + Knd) and 
where the Fourier components are developped in a Taylor series in power of the 
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small parameter e: 



l—p p—CO 

u(n,t) = YlJ2 e H l) (^,r n )A l (n,t), (6.4) 

l — — p p—1 
l—p p—OO 

p(n,t) = £ Y,^v r< £n,T n )A l {n,tl (6.5) 

l=- P p=\ 
l—p p=oo 

"(«.*) = E Y. eV € ) i^T n )A l (n,t). (6.6) 

The slow variables r„ and £„ are introduced via 

t„ = e(t -\ ) , f„ = e 2 n , (6.7) 

c 

where the velocity c will be determined later as a solvability condition. 

Using the identities 

A l (n+l,t)+A l (n-l,t) = 2cos(lKd)A l (n,t), (6.8) 

A l (n+l,t)-A l (n-l,t) = 2sin{lKd)A l (n,t), (6.9) 

the discret derivatives in (iO) can be written, with help of the derivative of a 



product (3.18), as 



-(u(n+l,t)-u{n-l,t)) = E E e^r/W^r^sm^d) 



l — —p p—1 

+ e(-)^-co S (Kld) 

C OT 

+ j[v P l) (rn + l,r) - 7j«(m - 1, r)] cos(Jfld) 

2 J £)2 (0 

+ ^(-) 2 ^r!- isin (^ w )M'K i ) ( 6 - 10 ) 



p— oo 



i(v(n+ l,t) - v(n - 1,*)) = E E £ p {4°K T ) Jsin (^ w ) 



i=—p p=i 

+ e(-)^cos(Kld) 

C OT 



~ [0»(m + l,r) - 0W(m - 1, r)] cos(Jfld) 

2 c ar^ 
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Subtituting ( |6.10 ), ( |6.11 ) in ( |6.3| ), we obtain 



l—p p—QO 



,(Q 



£ E e p {^ ) (™^)*sin(^d) + e(-)^cos(^d) 
i——p p—i 

£ 2 

+ -j [^°("» + l.r) - ??W(m - 1, t)] coa(Kld) 

2 j o2 CO 

2 c ar^ 

£=p p—oo 

- E E eP €W)^M) = o. 

i——p p—i 



(6.12) 



g P g £ P {e g^ I) + fflf^i) (TO> T)M i (n> t) 

l — — p p—1 

l—p p— cxo 

-E E eP ^ W)MW) = o, ( 6 - 13 ) 

l——p p—1 



l=p p—OO 



/i(0 



E E e P {^ ) (™^)*sin(KW) + e(-)^^cos(KW) 



l=—p P=l 



+ ir[<f>P{m + l,r) - ^W(m - 1, r)] cos(Kld) 



2 c ar 2 
" E E e ^ %, T) +'"^ ) Kr)}A'(n,t) 

l — —p p—1 

l=p p—QO l=p p—OO 

-i 1+ E E^WMVonE E c M I) ( m ' r M I (»'*)] = - ( 6 - 14 ) 
i——p p—i i——p p=i 

We proceed now to collect and solve different orders of e p and harmonics I 
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(order (p,l)) in ( 6.1 2| ) and (6.13) and (6.14). In the leading order (1,1) we have 

1=1 i=i 

Y^ r][ l) (m,T)ism(Kld)A l (n,t) ~~ ^ ip[ l) (m,T)A l (n,t) = 0, (6.15) 

i=-i i=-i 

i=i i=i 

^2 ini^{ l) (m,T)A l (n,t) - ^ <f>® (m, r) A 1 (n, t) = 0, (6.16) 

l=-i z=-i 

i=i i=i 

^ ^°(m,r)«sin(iffd)A J (n,*)- ^ iillT][ l) (m,T)A l (n,t) 



i=-i 



i=-i 



J2 i>['(m,T)A l (n,t) = 0. 



(6.17) 



i=-i 



Equations (3.15) ( |6.16 ) and (3.17) constitute a linear homogeneous system for 
rj[' ip[' and <j>\ ' . For I — equations ( |6.15 ) and (6.16) give 



V4 0) (™,t) =0, 
^ O) (m,r)=0, 



(6.18) 
(6.19) 



and (3.17) is satisfied. For I — 1 the determinant of the system for r][ , ip{ 
and <p\ is zero if 51 (K) verifies the dispersion relation 

sin Kd 



Cl(K) 



1 + sin 2 Kd 



(6.20) 



Unde r this c o nditio n we arrive at the following non-trivial solution of equations 

( FTHI) ( B (B) ' 



»?i (m,r) = J?(to,t), 

^ (iti,t) = i sm(Kd)i](m, r), 

(i) sin 2 Xd 

<Pi (w,t) = , 2 — T?(m,r), 

1 + sin A a 

where 77 is now the unknown function. 



(6.21) 
(6.22) 

(6.23) 
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At order (2, 1) we have the system 



1=2 



1 = 1 



£ vi\rn,r)i S in(Kld)A l (n,t) + £ (-) 



d drj[ '(m,T) 



l=-2 

1=2 



l=-l 



ch 



cos(Kld)A l (n,t) 



Y,4\™,T)A l (n,t) = 0. 



(6.24) 



l=-2 



1 = 2 



£ d ^ m ' T) A l (n,t)+ J2 Ml4 l \rn,T)A l (n,t) 



i=-i 

1 = 2 



ih 



l=-2 



J2<t>2 ] (m,T)A l (n,t) = 0, 



(6.25) 



(=-2 

1=2 



>-<- d aAQ 



V ^ !) (m,T)*sin(iadU'(n,<) + V (-)^-cos(K7dU'(n,i) 

Z — • Z — ^ r At 



l=-2 



l=-l 



1 = 2 



1 = 2 



1=1 Q (/)/ \ 

£ i W ' TJ ^^)-E ^°(m,r)}A'(n,*)- J] l/#W)A'(n,*) 






dr 



l=-2 

1 = 1 



l=-2 



£ ri[ l) (m,T)A l (n,t) x £ #(m,r)4'(n,f) = 0. 



i=-i 



i=-i 



For I — equation ( 6.25 ) gives using ( J6.18 ) 



4 0) (m,r)=0. 



(6.26) 



(6.27) 



Under ( 3.21 ) and ( 6.22 ), the equations ( 3.24 ) and ( |6.26| ) constitute a ho- 
mogeneous system for ip^ and the r derivative of r\ \ . The determinant is 
non-zero and consequently only the trivial solution exists, so that 

V4 0) (m,r) = 0, (6.28) 

(6.29) 



rf x (m, t) = 0. 



For I — 1 we have, using (6.25), an inhomogeneous linear system for rfe and 



(i) 



ip2 ' ■ The determinant of the associated homogeneous system is zero owing to 
the dispersion relation fl6.20| ). Therefore, the system will have a solution for the 
Fredholm solvability condition, which is satisfied for 



cos 3 (Kd) _ dVt 
(1 + sin 2 (Kd)) 2 ~ dK 

which determines c as the group velocity. 



(6.30) 



It is important to rematk here that, for the Taylor expansion (2.2), we must 
assume a nonvanishing group velocity c as indeed the first term in the expansion 
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of k(ui) is e/c. Hence we must avoid here the vicinity of the value Kd = ±ir/2 
for which c vanishes. In that vicinity, one should reconsider completely the 
problem. 

Then for all nonvanishing values of c we get 

r$\m,T) = g{m,r), (6.31) 

^ 1} (m, r) = ig(m, r) sm(Kd) + (j) ^K T ) cos (i«), (6.32) 

4 1} (m,r) = -n g (m, r) sm(ifrf) + t ^^ r) (sm(Kd) - tan ^\ ( 6 .33) 

where g(m,T) is an arbitrary function. Furthermore for I = 2 we obtain ^ 
and ip 2 an d ^2 as follows: 

r] 2 2 \m,T) =a(Kd)r) 2 (m,T), (6.34) 

V>2 2) (to,t) = isin(2Kd)a(Kd)r 1 2 (m,T), (6.35) 

4 2) (m,r) = -20sin(2 J ft:d)a( J ft:d)r7 2 (m,r), (6.36) 

where a(Kd) is defined as 

a(Kd) = ^S -. (6.37) 

20 [1 + sin 2 (2Kd)] + sm(2Kd) 

The next order (3, I) gives the system 

1=3 1=2 , „ (l) 

J2 V3\m,T)ism(Kld)A\n,t)+ J^ (~)^^ cos(Kld)A l (n,t) 

l=-3 l=-2 



1=1 

2 

!=1 i w «2„(0 



^ -[^(m + l.r) -r)[ l) (m- l,T)]cos(Kld)A l (n,t) 



S ^) 2 ^^in(^d)^(n,i) 



z=-i 

i=3 



24 ("»,r)A'(n,t) = 0. (6.38) 



2=-3 



£ ^ ; W ' T) AW)+E i«^)(m,r)M'(n,«) 



<9t 

i=-2 i=-3 

1=3 



^ 4' ) (m,r)}A'(n,i) = 0, (6.39) 



i=-3 
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/=3 l=2 r\ 8rh^( \ 

J2 4° (m,T)ism{Kld)A l (n,t) + ^(-) feim ' T) cos(Kld)A l (n,t) 

l=-3 l=~2 

1 = 1 

+ J2 ^i ) {m+l,T)~(f ) ( p(m-l,T)}cos(Kld)A l (n,t) 
i=-i 

Z=2 q (J)/ \ '=3 

- E %7' T) AW)-]T tfM^(m,r)M'(n,t) 

i = -2 (=-3 

1=3 1=1 1=2 



]T Vf (m,r)A l (n,t)- J2 v[ l \m,T)A l (n,t) x ]T $\m,T)A l (n,t) 

= -3 l=-l l=-2 

1=2 1=1 

Y^ v 2 l) ( m > t ) a ' («> i ) x Yl ^i° ( m ' t ) a ' ( n ' *) = °- ( 6 - 4 °) 



l = -2 l=-l 



The order (3, 0) determines rfo , ty$ and </> 3 via an inhomogeneous system 
of equations. They read 

^( m , T ).-^M|,( mi r)r, (6.41) 



^Kr)^(-)^I#.V!P. «,42) 



c rf + c 
4 0) (™,t)=0. (6.43) 

The next order (3, 1) allows us to find the nonlinear evolution of r\. It is an 

inhomogeneous linear system of equations for r\§ , V4 an d 4>s of determinant 
zero. Therefore, we will have a solution if the Fredholm solvability condition 
holds. This condition gives the nonlinear evolution of 77 in which the term in 
g(m, t) coming from rfe , V4 anQl $2 cancels out. The equation for 77(771, r) = 
r/ m reads finally 

f)2 

-if3(Vm+i - Vm-i) + a ^r~ ~ 7l'7m| 2 »7 m = , (6.44) 
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with the following definitions (we set Kd = K for simplicity) 

1 cos 3 K sin K 
~2 1 + sin 2 if ' 

1 (cos 2 if -l)(cos 2 if + 6)(cos 2 if-2) 2 



a = — 



2 cos 4 K 

1 cos 2 K -2 



2 cos 2 if -4 

8 cos 6 K + 6 cos 5 if - 20 cos 4 if - 11 cos 3 if + 13 cos 2 if - 4 

V 

4 cos 3 if + 3 cos 2 if - cos if + 1 

Note that inn the Brillouin zone if G [— 7r, +7r] , the coefficient a has a singularity 
in if = ±7r/2 where the group velocity vanishes, which is a forbidden region. 



7 Three wave resonant interaction 

To further illustrate the method, we consider here the Takeno discrete model 
for interaction of excitons (or vibrons) with the phonons in a lattice of coupled 
harmonic oscillators (via a Fr61ich-like hamiltonian) Jl7p . The model results 
from the Hamiltonian 

H = Hph + W cx + Wint > (7.1) 

H p h = - ^ Mill + S[u n+1 - u n } 2 



Wcx = -Z ^2 m \in + u oQn] + »[?n+l ~ Qn] 



Hint = X X! ^[ W "+l ~~ u n-l]g£ ( 7 ' 2 ) 



After the rescaling 



q n = —V2mMq' n , u n = — u^ (7.3) 



and forgetting the primes, the equations of motion are 

u n - Vl\[u n+ i - 2u n + u n -i] = ql + i - Q 2 _i (7.4) 

Qn - <*>i[g n +l - 2 Qn + Qn-l] + ^o<7™ = -Qn[u n +1 ~ «n-l] 



with the following definitions 



n\ = — , u; 2 = - 
1 M ' * m 
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7.1 Resonant wave interaction 

We consider now the situation of 3-wave scattering, that is the situation where 
the exciton wave contains two components (frequencies v\ and 1/2) which interact 
with the phonon wave (frequency Q) according to the Brillouin selection rule 

v\ — V2 = £1 , fci — hi = K . (7-5) 

For slowly varying envelopes we set 

q n {t)=e[a { f ) e i ^ + a 2 1) e i(> *] (7.6) 

+ afe^+V + afe^-V] + c.c. + 0(e 3 ) , 
u n (t) =ebV>e** (7.7) 

+ e 2 [bf>e 2 ^ + bfe™* + bfe 2i * + bfe* 9 ^ 
+ b f e m+V + h f e HH~^] + c.c. + Q(e 3 ) , 

where the amplitudes af and of depend on the slow variables (m, r) and with 
the definitions 

t = et , {£„ = en e = — } -> m , 

B\ = kin — v\t , 82 = k2n — V2t , <j) — Kn — fit , (7-8) 



The choice of the above harmonics in the order e 2 results from the remark 
that quadratic terms (like q n u n ) for 3 waves induce (with the complex conju- 
gates) waves with phases 29\, 28 2 , 2(f), 81+82, @i + 4>, $2 + 4>> Q\ — $2, #i — 4>, 
8 2 — <t>- But the selection rules give 8 2 + <f> = 9\, 6\ — 8 2 = 4>, 81 — <f> = 8 2l which 
are already considered at first order and hence need not to be included in the 
second order. 

The above change of variables has to be applied to functions 

for which, using the tools developped in sec. now for N = e _1 , we readily 
obtain at order e 

Cp n = [-w 2 Vm - e2iu>d T iP m + ■ ■ -] e l{kn -" *>, 

ip n+1 -<p n -i = [2isink i\) m + e cos k(if) m+ i - V'm-i) H ] e l (fc™-^) ; 

ip n+1 - 2tp n + ip n _i = [2(cosk-l)ip m + iesmk(tp m +i - V'm-i) H yV™-"*) , 



2S 



7.2 The limit equation 



All the above machinery is applied now to the system (7.4) which gives at order 
e : 

e i<t> . [-n 2 - nl2{cos K - l)]b ( ^ = 0, 



tt>i 



2l„(l) 



-v{ - w{ 2(cos fci - 1) + WoK = 0, 



e je2 : [-i/f - w 2 2(cos fe 2 - 1) + w^a^ = 0. 

These imply the dispersion relations 

■>2 or>2/ 



fl z =20f(l- cos iT) , 



v\ = u>q + 2w 2 (l - cosfci) , 
and hence the 3 group velocities 



Vo 



u>, 







•w 2 (l — cos/c 2 ) 



fi? . 



n 



sin if 



Vl 



UJ 



V\ 



- sin k\ 



v 2 



— sin k 2 



V 2 



The next order, e 2 gives in turn 



i0 2 



2^ 1} -fi^siniT^ 



(i) 

m+1 

(i) 



^(1) 



e 

„2i0 



„2i6>. 



„2id 



o i(Sl+0) 



oWi-4>) 



i l,m 



] = ^sinif a^ a 2 , 
..J = -2isini<: a^&f', 



(7.9) 



(7.10) 



-2iv>ia\ — 0Jii$m.ki[a\ „ l+1 

-2iv 2 a 2 — ujlisiiik 2 [a 2 ' m+1 — ai^-i] = 218111^ a 1 i; 6 1 i; , 

-4iy 2 -2fi 2 (cos2fci- 1)] b (2) = iisinh cosfci [a^] 2 , 

-4zy 2 - 2cj 2 (cos2fci - 1) H-Wq] a^ 2) = 0, 

-4i/| - 20 2 (cos2fc 2 - 1)] 4 2) = 4isinfc 2 cosfc 2 [aj 2) ] 2 , 

-4^ 2 - 2w 2 (cos2fc 2 - 1) + ujI] a 2 2) = 0, 

-4fl 2 - 2Q 2 (cos2if- 1)] bf ] = 0, 

-40 2 - 2cj 2 (cos 2K-1) + Wq] 4 2) = °> 

-4(^i + ^ 2 ) 2 - 2ft 2 (cos(fci + fc 2 ) - l)]6f ) = 4isin(fc 1 +fc 2 )a 1 1) 4 1) , 

-A{ Vl + v 2 ) 2 - 2w 2 (cos(fci + k 2 ) - 1) + w 2 ] 4 2) = 0, 

-4fa + 4>f - 20 2 (cos(fci + K) - 1)] ^ 2) = 0, 

-4(z/i+<£) 2 - 2w 2 (cos(fci+^) - 1) + Lu 2 ]a ( 5 2) ■- 

-A{v 2 - 0) 2 - 20 2 (cos(fc 2 -K)- 1)] bf ) = 

-4(iA2 - </>) 2 - 2w 2 (cos(fc 2 - K 



2ismK a^ ] b^\ 



i)+c 2 ]4 2) 



2isinK a 2 1] b { i ] . 



The above first 3 equations provide the evolutions of the envelopes which we 
scale as 



',<D 



X = sinK b\ , a\= svnK a\ 



(i) 



• V C 1 ) 

a 2 = sin K a 2 



2') 



and consequently which obey 

1 2 

X + v-[X m+1 - X m _i] =—q a i^2 , 

61 + vi-[ai l7n +i - aim-i] = — a 2 X , (7-11) 

a 2 + v 2 -[a2,m+i ~ a 2 . m _i] = a 1 X . 

I v 2 

The remaining 12 equations give simply the coefficients of the second har- 
monics in terms of those of the first, precisely: 

(2) (2) ,(2) (2) (2) ,(2) ,(2) n 

<A = a 2 = b 3 = a 3 = a 4 = b 5 = b 6 = ° > 

,(2) _ . I/? + ^nf(cOs2fci-l)] (!) 2 

sin k\ cos k\ 
,(2) _ . v\ - |»i(cos2fc 2 -1) (2) 



b\ A > = 



sin k 2 cos k 2 

\2 1 lo2 



(2) _ . Ql + ^2) + |nj(cOS(fcl + feg) - 1) (1) (i) 

sm(fci + fc 2 ) 

(2) _ , 2(1/! + 4>f + C^COSJfci + X) - 1) + a;g (!)(!) 

5 smK 1 x 

(2) _ . 2(^ 2 - 0) 2 + ^ 2 (cOs(fc 2 -K)-l)+ul (!)-(!) 

6 smK 2 x 



Hence eqs. fl7.1l| ) constitute a closed system of equations for the first order 
variations of the envelopes. It is the discrete analogue of the continuous 3- wave 
resonant interaction system. 



A Discrete change of variable 

For a given function x of the discrete variable n, we are intrested in expressing 
all-order differences in a large grid indexed by the variable m (taking points sep- 
arated by a given odd interger p) in terms of the original hierarchy of differences 
([03). The following notations will be used throughout 



V l x n = J2(-) k C\ x n+e - 2k = a#> , (A.l) 

fc=0 

1 

Ap£„ = $^(-)* C t X n+p(l-2k) = x<& . (A.2) 



fc=0 



To express the A-differences in terms of the V-differences, we first write exact 
Taylor like series as the following identities (valid for any given n and any point 
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function x n ) 

P ■ ■ 1 P 1 

x n+p = / j a p x n ' 7y~ x n -i x n—p — / j \~ ) a p x n ' \) 7} X n i (,-"-. oj 

3=0 ' j=0 

with the following definitions for the tiWe-derivatives 



p 
n+ P + £(-)' C e p+1 x n+p+1 _ 2e + (-) p+1 x n ^ p . (A.4) 

=1 



x (p) = x -+- W-V ^ . . *_ -4. r_^+! 



For instance we may write 



r _,_- -r + -r^ 4- -t< 2 ) 4- -r< 3 ) + -t< 4 > + -x^ + if( 5 ) 

(1) (2) ~(2) 

^n+2 =X n + —X n + —X n + "I„ 

• T n+1 =3; ri + 7y~ x n • 7} x n 

X n —1 =x n ~ 7y~ x n • 7} x n 

X o =x - - X W + -x^ - -x^ 

■r — rp _ _ 7 .(1) i t(2) 7.(3) _i ™(3) 

r A -r - -rW 4- -t< 2 ) - -r< 3 ) 4- -rW - -r< 4 ) 

x * -a; - -xW 4 -x( 2 ) - -x<3) + -i<") - -a;( 5 ) 4 -x^ 



which are useful to play with in order to understand in a concrete way the 
relations that follow. 

The a J p are the coefficients to be computed but only some of them are needed, 
as indeed we consider differences A p x n - For the definitions 

a{ = 2af q + + \, # = 2o£, (A.5) 

the first and second derivatives read (p = 2q 4- 1) 

<j 

Xn+p X n — p > 0;^ V X n , ^A.DJ 

r 
X n+ 2r - 2x n + X n -2r = "Y] &r V 2£ X„ , (A. 7) 
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and we prove hereafter the expressions ( p.8| ) and (3.£) for the coefficients a q and 
ft 0: 



t (2g+l)(g + l)! e 2r(r + £-l)\ 

<X„ = 7TT , P r = 



(A.. 



9 (q - £)l(2£ + 1)\ ' " r (r-£)\(2£)\ 
Note that in ( |3.9[ ) we simply have set J q = ft for r = 2q + 1 . 

As the point n is fixed (arbitrary) lighter notations can be used, namely 



n - T (2«+i) 
^q — ^n 

H -x {2r) 



Aq — X n +(2q+l) ~ %n— (2g+l) ; 
<_>0 = 2-n i L>r = ^n+2r T 2^n— 2r i 



and the definition (A.l) can now be written 

q r 

G 1 = EH 9 "" C lq+1 A * ' ff r = EH r ~ fe C5T* S fc 



fc=0 



fc=0 



while (A. 6) and (A. 7) become respectively (note that /3" = 2) 



A q =J2a e q G e , S r =J2tfH r . 



(A.9) 



(A.10) 



(A.ll) 



£=0 



Note that neither in A q nor in S r appear the tilde-differences (A.4). 

The first step is the computation of the coefficients a q . By replacing ( A.fO) 
in (A.ff ), we arrive at the following equivalent relations 



^=e«:eh"c^ 



£=0 fc=0 

9 



G 9 = E(-)^ fec U £«to 



(A.12) 



(A.13) 



fc=0 



fco 



which are identities valid for any function x n , that is for any choice of the 
sequences {Ak} and {Ge}. Consequently the coefficients of each Ak and Gk 
identically vanish, which furnishes the equivalent recursion relations 



V{^, 9 }, £< q , j(-) w ^aj = 



1 



aj = l 



(A.14) 
(A.15) 



fc=^ 



The rest of the proof is the check that the expression of oc„ given in (A.8) 



does solve the above recursion relations. This is easily done with help of the 
following identity 



m=0 



, (s + n + m — 1)! 
m!(n — to)!(s + to)! 







(A.16) 
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which can be obtained by differientiating n times a conveniently chosen poly- 
nomial of degree n — 1 in the variable s, namely 

n 

Df(s) = f(s + 1) - /(,) , D n f(s) = J2 (-) n - m C™ f(s + m) , (A.17) 
(s + n-2)! 



/(* 



(s-1)! 



D n f(s) = 



(A.18) 



The same procedure is then applied to get the coefficients fi^ as given in (A.8). 

Aknowledgements: It is a pleasure to aknowledge many helpful discus- 
sions with G. Mennessier whose contribution has been essential in many tech- 
nical points. 
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